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Context. We use Hubble Space Telescope photometry of six rich, compact star clusters in the Large Magellanic Cloud (LMC), with 
ages ranging from 0.01 to 1.0 Gyr, to derive the clusters' stellar mass functions (MPs) at their half-mass radii. 

Aims. The LMC is an ideal environment to study stellar MPs, because it contains a large population of compact clusters at different 
evolutionary stages. We aim to obtain constraints on the initial MPs (IMPs) of our sample clusters on the basis of their present-day 
MPs, combined with our understanding of their dynamical and photometric evolution. 

Methods. We derive the clusters' present-day MPs below LOM© using deep observations with the Space Telescope Imaging 
Spectrograph and updated stellar population synthesis models. 

Results. Since the relaxation timescales of low-mass stars are very long, dynamical evolution will not have affected the MPs below 
I.OMq significantly, so that - within the uncertainties - the derived MPs are consistent with the solar-neighbourhood IMP, at least for 
the younger clusters. 

Conclusions. The IMP in the low-density, low-metallicity environment of the LMC disk is not significantly different from that in the 
solar neighbourhood. 
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1. Introduction 

The shape of the stellar initial mass function (IMF) is a very 
important unresolved issue in modern astrophysics, because it 
plays a crucial role in many of the remaining 'big questions.' 
The IMF is usually assumed to be universal, and best approx- 
imated by either a power-law (e.g., Kroupa 2001) or a lognor- 
mal distribution (Chabrier 2003; Andersen et al. 2008). Kroupa 
(2001) studied the Galactic-field IMF down to 0.0 IMq and de- 
rived a three-part power-law function. Chiosi et al. (2007) ob- 
tained a similar mass function down to O.IMq based on their 
analysis of three clusters in the Small Magellanic Cloud (SMC), 
and concluded that the IMF in SMC clusters is in agreement with 
the 'standard' Kroupa (2001) solar-neighbourhood IMF. Da Rio 
et al. (2009) recently studied the stellar association LH 95 in 
the Large Magellanic Cloud (LMC), focusing on its pre-main- 
sequence (PMS) stars. Their results showed that there are no sig- 
nificant differences between the IMFs of the entire LH 95 region 
and those of three individual subclusters: they all follow a mul- 
tiple power-law distribution. Andersen et al. (2008) studied the 
low-mass stellar mass distributions of seven star-forming regions 
and concluded that the composite IMF is consistent with a log- 
normal distribution. Paresce et al. (2000) obtained a similar re- 
sult after analysing the mass functions (MFs) of a dozen Galactic 
globular clusters (GCs) for stellar masses below I.QMq. Thus, 
the form of the universal IMF appears to be best approximated 
by both Kroupa (2001) broken power-law and lognormal distri- 



butions (Chabrier 2003; Andersen et al. 2008; Covey et al. 2008; 
Liu et al. 2009; Oliveira et al. 2009). To distinguish between ei- 
ther shape, one would need to probe down to the stellar/brown- 
dwarf transition region, which still poses a significant observa- 
tional challenge, particularly in (even the nearest) extragalactic 
environments. 

The evolution of the stellar MF is also important more gen- 
erally, because the IMFs of many clusters and galaxies cannot 
be observed directly. However, we can derive their IMFs based 
on the present-day MF (PDMF) if - at least - we understand 
its evolution in detail. Star clusters, both open clusters and GCs, 
provide ideal objects to tackle many astronomical problems, be- 
cause all of their member stars have approximately the same age 
and metallicity, and are located roughly at the same distance. 
Although the GCs in the Milky Way are relatively nearby and 
their members can be observed easily, they are not well suited to 
study the evolution of the MF, because Galactic GCs are all old 
(with ages f > 10 Gyr). They can therefore only provide evolu- 
tionary information on long timescales. Galactic open clusters, 
on the other hand, are only effective tracers of MF evolution on 
short timescales, while they also tend to be affected quite sig- 
nificantly by small-number statistics. Ideally, therefore, we need 
rich massive clusters covering a large age range to make signif- 
icant progress on this important problem. This makes the LMC 
an ideal laboratory, because it contains a large population of rich 
star clusters with masses similar to Galactic GCs and covering 
ages from 0.001 to 10 Gyr (e.g., Beaulieu et al. 1999; Elson et 
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al. 1999). This implies that we can study the MF at almost all 
evolutionary stages using the rich, compact clusters in the LMC. 
In the past, it proved impossible to resolve individual stars in 
dense star clusters at the distance of the LMC (~ 50 kpc), but 
the unprecedented, high spatial resolution of the Hubble Space 
Telescope (HST) facilitates such studies today. 

This is what we set out to do in this paper In Section 2, we 
briefly describe our HST observations and give a basic overview 
of the data-reduction procedures (see detailed steps in Liu et al. 
2009). We present and discuss our main results in Sections 3 and 
4, respectively, and provide a summary in Section 5. 



2. Observations and data reduction 

2.1. Observations and previous work 

As part of HST programme GO-7307 we observed a carefully 
selected cluster sample in the LMC, including six compact clus- 
ters in three pairs (Pair I: NGC 1805 and NGC 1818, Pair II: 
NGC 1831 and NGC 1868, and Pair III: NGC 2209 and Hodge 
14): see Table 1 for their fundamental parameters. Our six sam- 
ple clusters have ages of 10^ - 10^ yr, with the additional con- 
straint that the two clusters in each pair have similar ages, metal- 
licities, total mass, and distance from the LMC centre, yet dif- 
ferent structural parameters (Beaulieu et al. 1999). 

Much work has been done already based on these high- 
quality imaging observations (e.g., Elson et al. 1998, 1999; 
Johnson et al. 2001; Santiago et al. 2001; de Grijs et al. 
2002a,b,c; Kerber et al. 2006; Liu et al. 2009). de Grijs et al. 
(2002a,b) obtained the MFs of the two youngest clusters in our 
sample for stellar masses above 1 .QMq, and concluded that they 
are largely similar to the Salpeter (1955) IMF. Liu et al. (2009) 
analysed the MF of NGC 1818 below 1.0 Mq and combined 
their low-mass MF with the higher-mass results of de Grijs et al. 
(2002b). They found that the IMF of NGC 1818 could be wefl 
approximated by both a Kroupa (2001)-type broken power-law 
function and a lognormal distribution. 

Our high-quality imaging observations were obtained with 
both the Wide-Field and Planetary Camera-2 (WFPC2) and the 
Space Telescope Imaging Spectrograph (STIS). WFPC2 is com- 
posed of four chips (each containing 800 x 800 pixels), one 
Planetary Camera (PC) and three Wide-Field (WF) aiTays (WF2, 
WF3, and WF4). The pixel size of the PC chip is 0.0455 arcsec 
(with a field of view of ~ 34 x 34 arcsec^) and that of each WF 
chip is 0.097 arcsec (with fields of view of ~ 150 x 150 arcsec^). 
The STIS pixel size is 0.0507 arcsec and the corresponding field 
of view is about 28 x 52 arcsec^. 

WFPC2 exposures through the F555W and F814W filters 
(roughly coiTesponding to the Johnson-Cousins V and / bands, 
respectively) were obtained with the PC chip centred on both 
the clusters' half-mass radii (with a total exposure time of 2500 
s in both filters; see for more details Santiago et al. 2001; de 
Grijs et al. 2002a). Deep STIS exposures in ACCUM imag- 
ing mode through the F28x50LP long-pass filter (central wave- 
length ^ 7230A) were also obtained, centred on the clusters' 
half-mass radii (with total exposure times of 2950 s for Pairs I 
and III, and 2890 s for Pair II; see Elson et al. 1999). To ob- 
tain 'clean' MFs, we must subtract the background stellar con- 
tribution. We therefore also obtained deep WFPC2 images from 
the HST Data Archive of the general LMC background through 
the F555W and F814W filters, with exposure times of 7800 and 
5200 s, respectively (see Castro et al. 2001; Santiago et al. 2001; 
de Grijs et al. 2002a; Liu et al. 2009). 
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Fig. 1. Colour-magnitude diagrams of our LMC clusters and the 
best-fitting Padova isochrones (Girardi et al. 2000). The hori- 
zontal dashed and solid lines in each panel represent the upper 
and lower magnitude limits, respectively, of the parameter space 
covered by our STIS observations. 



2.2. Data reduction and pliotometry 

We use the same method for data reduction and photometry as 
in Liu et al. (2009). Aperture photometrjQ was performed on our 
WFPC2 and STIS images using the iraf/apphotEI package, and 
2-pixel apertures were adopted since this produced the smallest 
photometric errors and the tightest cluster main sequences. 

We used the relations of Whitmore et al. (1999) to cor- 
rect the resulting photometry for the time-dependent charge- 
transfer efficiency effects and iraf/stsda^ to rectify the geo- 
metric distortions of the WFPC2 chips. We subsequently ap- 
plied aperture corrections based on the model PSFs generated by 
Tiny Tim (Krist & Hook 2001). Finally, we used the relations of 
Holtzman et al. (1995) to convert the aperture-corrected F555W 
and F814W magnitudes to the Johnson-Cousins V and / pass- 
bands. Figures[T]and|2]show, respectively, the resulting WFPC2- 
based colour-magnitude diagrams (CMDs) of all sample clusters 
and the spatial distributions of the stars in our LMC cluster fields 
as observed with both WFPC2 and STIS (dots and rectangular 
outlines, respectively). 



' See Liu et al. (2009) for a discussion of the pros and cons of using 
aperture photometry versus point-spread function fitting. 

' The Image Reduction and Analysis Facility (iraf) is distributed by 
the National Optical Astronomy Observatories, which is operated by 
the Association of Universities for Research in Astronomy, Inc., under 
cooperative agreement with the US National Science Foundation. 

^ STSDAS, the Space Telescope Science Data Analysis System, con- 
tains tasks complementary to the existing iraf tasks. Version 3.1 was 
adopted for the data reduction performed in this paper. 
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Cluster logCage)" E{B-Vf {m-M)l \og{M,,IM^y D^^^ 

(pc) (°) 

NGC 1805 7.65±0.05 0.04 18.59 3.52+0.13 1.33±0.06 3.86-4.00 

NGC 1818 7.65±0.05 0.03 18.58 4.13!;; 2.45±0.09 3.47-3.61 



[yr] 


(mag) 


(mag) 


7.65±0.05 


0.04 


18.59 


7.65±0.05 


0.03 


18.58 


8.75±0.05 


0.00 


18.58 


9.00±0.05 


0.02 


18.55 


9.10±0.05 


0.07 


18.39 


9.30±0.05 


0.04 


18.49 



NGC 1831 8.75±0.05 0.00 18.58 4.81±0.13 4.44±0.14 4.82-4.85 

NGC 1868 9.00±0.05 0.02 18.55 4.53±0.10 1.62±0.05 5.57-5.47 

NGC 2209 9.10±0.05 0.07 18.39 5.03!!!-^'' 5.43±0.33 5.48-5.43 

Hodge 14 9.30±0.05 0.04 18.49 4.33!^^^ 1.80±0.14 4.19^.37 

Notes: Dlmc is the distance from the LMC centre; the two values indicate Dlmc to the optical, geometrical centre (Bica et al. 1996) and the 
dynamical, rotation centre (Westerlund 1990), respectively. References: " this paper (the age uncertainties are driven by the discreteness of the 
Padova isochrones), * Castro et al. (2001), Mackey & Gilmore (2003), Meurer, Cacciari, & Freeman (1990). 
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Fig. 2. Spatial distribution of the stars in our LMC clusters. The 
dots represent the stars detected in the WFPC2 observations, 
while the solid rectangular shapes indicate the areas covered by 
our STIS observations, which are used for further analysis in this 
paper. 



2.3. Completeness corrections and background subtraction 

One of the most difficult problems we faced to derive clean MPs 
involved correcting for sampling incompleteness, which is usu- 
ally a function of position in a cluster We used the same method 
as in Liu et al. (2009), in essence a slightly modified version of 
the approach used in de Grijs et al. (2002a), who computed these 
connections in circular annuli around the cluster centres. We 
computed the completeness corrections for the entire STIS chip, 
because the STIS observations were centred on the low-density 
half-mass radius and the effects of sampling incompleteness are 
constant across our STIS field within the observational uncer- 
tainties (e.g., de Grijs et al. 2002a; Liu et al. 2009). We used 
the same method to compute the equivalent completeness cor- 
rections for the background field. We added an area-dependent 
number of artificial sources of Gaussian shape to each STIS ex- 
posure with input magnitudes between 16.0 and 30.0 mag, in 
steps of 0.5 mag. We then adopted the same photometric anal- 



ysis method for the fields including both the 'real' cluster (and 
background) stars and the artificial stars, to asses how many ar- 
tificial stars could be detected after correction for blends and 
superpositions. For the analysis in this paper, we only consider 
magnitude ranges that are > 50% complete. 

To obtain a clean MP, we must subtract the contamination by 
the background field, because its stellar mass distribution is gen- 
erally different from that of the clusters (cf. Castro et al. 2001). 
We do not have a background field in the STIS F28x50LP pass- 
band. Instead, we used a general WFPC2/F814W background- 
field observation. However, because our STIS observations are 
much deeper than the WFPC2 data, we could not directly 
subtract the full background effects from the observations. 
Gouliermis et al. (2006a) suggested that the stellar mass distri- 
bution in the LMC disk follows a broken power-law distribution. 
We therefore adopted a power law (with slope F = 1.87, where 
the IMF, ^(m») oc trf) to approximate and extrapolate the stellar 
mass distribution in the general LMC background down to 0. 1 
Mq (see Liu et al. 2009, their fig. 5). 

In Figs.[3]and|4]we show the photometric completeness frac- 
tions and the resulting luminosity functions (LFs) for all of our 
sample clusters, respectively. The LFs have been fully complete- 
ness corrected and background subtracted. We only include mag- 
nitude bins for which the observational completeness fractions 
are greater than 50%. 

3. Analysis 

3. 1 . Age, metallicity, and evolutionary models 

Santiago et al. (2001), de Grijs et al. (2002b,c), and Kerber et 
al. (2006) studied the MFs of our sample clusters above I.OMq, 
based on WFPC2 data. Liu et al. (2009) studied the MF of NGC 
1818 below 1.0 Mq based on STIS data. Because the STIS data 
are much deeper than the WFPC2 measurements (e.g., de Grijs 
et al. 2002a), here we study the MFs of all six LMC clusters from 
HST programme GO-7307 based on the STIS observations. 

The ages of our six sample clusters cover a large range, from 
10^ - 10^ yr. Most of the stars below l.OM© of the clusters in 
the youngest of pair, NGC 1805 and NGC 1818, are still on 
the PMS (Liu et al. 2009). The evolution of PMS stars is still 
uncertain (Baraffe et al. 1997, 1998). White et al. (1999) con- 
cluded that the models of Baraffe et al. (1997, 1998) resulted 
in the most consistent ages and masses, on the basis of a com- 
parison of six PMS evolutionary model sets (Park et al. 2000). 
Because our STIS observations were obtained in a single pass- 
band only, we cannot derive the PMS ages and metallicities on 
the basis of CMD analysis. However, the cluster stars from our 
WFPC2 and STIS observations occupy a common locus on the 
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Fig. 3. Completeness ratios of all sample clusters. The dashed 
lines represent the 50% completeness limits for the individual 
clusters. 
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Fig. 4. Completeness-corrected, background-subtracted LFs of 
all sample clusters, for completeness fractions > 50%. 



CMD, so we can obtain the basic cluster parameters by adopt- 
ing Padova isochrones (Girardi et al. 2000) to fit the faint end 
of the main sequence, and the Baraffe et al. (1997, 1998) mod- 
els for the sequence of low-mass stars, i.e., by extrapolating the 
WFPC2 observations to fainter luminosities (see, for details, Liu 
et al. 2009). 

The other sample clusters are old enough for all low-mass 
stars to have evolved onto the main sequence, so that we can get 



their ages and metallicities using Padova isochrones for main- 
sequence fitting. The best fits and the resulting basic cluster pa- 
rameters are shown in Fig.[T] We used the Padova isochrones to 
fit the main-sequence ridge line in the CMD with diff'erent metal- 
licities and ages for each cluster Next, we compared the quality 
of all fits and adopted the best-fitting metallicity and age for the 
cluster of interest. The small spread of the main-sequence data 
points will cause a small uncertainty in the age determination 
for each cluster, which we characterise by using the discrete- 
ness of the Padova isochrones (in steps of 0.05 dex in age) as a 
proxy for the age uncertainties. In fact, for each cluster the near- 
est isochrones to the best-fitting model provide markedly worse 
fits to the main-sequence turn-off location, so that we are confi- 
dent that our fits are robust. 

In Liu et al. (2009) we added the F28x50LP filter to the 
Baraffe et al. (1998) model suite. These models cover a range 
of metallicities, enabling us to choose the most relevant mass- 
luminosity relation for conversion of the F28x50LP magnitudes 
to individual stellar masses. 

A significant body of work exists in the literature to support 
our choice of metallicity for the individual clusters. For NGC 
1805, Johnson et al. (2001) used HST CMDs to derive a mean 
[Fe/H] ~ (solar metallicity), although MeHani et al. (1994) had 
argued previously that the cluster's most appropriate metallicity 
was Z = 0.008 (where Z© - 0.020), i.e., the average metallicity 
of the young LMC field population. On this basis, we adopted 
Z - 0.008 for this cluster, corresponding to [Fe/H] ^ -0.37 
(assuming a one-to-one correlation between metallicity and iron 
abundance). Using solar metallicity instead, we derive an age for 
NGC 1805 of log(r/yr) ^ 7.50 + 0.05 using Padova isochrones. 
While the slopes of the resulting MFs based on either metal- 
licity are similar (see Fig. |5] top), the calibration of the solar- 
metallicity MFs would shift to higher masses. 

NGC 1818 has been the subject of a large number of stud- 
ies aimed at determining its metallicity. Although Johnson et 
al. (2001) obtained solar metallicity from HST CMD analy- 
sis, most other modern (predominantly spectroscopic) deter- 
minations centre around either [Fe/H] - -0.37 ± 0,03 (e.g., 
Jasniewicz & Thevenin 1994; Bonatto et al. 1995) or [Fe/H] 
~ -0.8, roughly corresponding to Z = 0.003 (e.g., Meliani 
et al. 1994; Will et al. 1995; Oliva & OrigHa 1998). The MF 
slopes based on a metallicity of Z = 0.008 are steeper than 
for Z - 0.004 (the lowest metallicity isochrone available for 
the cluster's young ag^ (see Fig. \5\ bottom), although they 
are both in good agreement with the Kroupa (2001) IMF slope 
(F = 0.3 + 0.5). 

The most recent metallicity estimates for NGC 1831 con- 
verge to [Fe/H] ~ -0.35 (Bonatto et al. 1995, based on UV spec- 
troscopy), while Vallenari et al. (1992) similarly suggested a best 
estimate of Z = 0.008 based on their analysis of the literature on 
this cluster at the time of their publication, but see Mateo et al. 
(1987) and Olszewski et al. (1988, 1991) for close-to-solar abun- 
dance estimates. Given the current observational status for this 
cluster, we adopted Z - 0.008 as a compromise. 

Both NGC 1868 and Hodge 14 are somewhat more metal- 
poor than the younger clusters in our sample. For NGC 1868, 
Bica et al. (1986) reported [Fe/H] = -0.6 + 0.35, consistent 
with Olszewski et al.'s (1991) spectroscopic metallicity deter- 



Although the Padova models include isochrones for Z = 0.001, 
they are only provided for ages in excess of log(//yr) = 7.80. Use of the 
Z = 0.004 isochrone results in a cluster age of log(?/yr) = 7.65 ± 0. 10. 
For the low-mass PMS stars, we use Z = 0.003 and log(f/yr) = 7.25 
based on the Baraffe et al. (1997, 1998) models. 
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mination, [Fe/H] = -0.50. Similarly, Jensen et al. (1988) and 
Olszewski et al. (1991) used spectroscopy of Hodge 14 to obtain 
[Fe/H] = -0.66 + 0.2. For both clusters we adopted Z = 0.004, 
corresponding to [Fe/H] - -0.68. 

Finally, NGC 2209 is the lowest-metallicity cluster in our 
sample, with [Fe/H] = -0.9 + 0.3 (Bica et al. 1986; see also 
Chiosi et al. 1986; Dottori et al. 1987; Frogel et al. 1990). Based 
on a careful analysis of the goodness-of-fit parameters for the 
range of metallicities provides by the Padova isochrones, we 
adopted Z = 0.004 ([Fe/H] = -0.68) for NGC 2209; its CMD is 
much more poorly approximated for [Fe/H] - -0.9. 

Metallicity does not play an important role in our analysis 
of the MF slopes. The models of Baraffe et al. (1997, 1998) for 
different metallicities yield similar main sequences at low mass. 

3.2. Mass functions 

We adopted the same method as Liu et al. (2009) to derive the 
MFs of the clusters of pair L Although the age of NGC 1805 and 
NGC 1818 is log(f/yr) ^ 7.65, many low-mass stars are still on 
the PMS. The ages of the PMS stars in both clusters are about 
log(f/yr) ^ 7.25 + 0.40 (Liu et al. 2009). All stars in the clusters 
of pairs II and III have already evolved onto the main sequence. 

Several studies explored the effects of mass segregation in 
these clusters (based on the same WFPC2 observations used 
here) by dividing the full field of view into a number of smaller 
areas at a range of distances from the cluster centres (e.g., 
Santiago et al. 2001; de Grijs et al. 2002a,b,c; Kerber et al. 
2006). However, our STIS observations were taken at the half- 
mass radii of our sample clusters and the STIS field is much 
smaller than that of the combined set of WFPC2 observations, 
so we limited our analysis to both the entire STIS region and 
areas at two different radii. Figures |5] |6] and |7] show the cluster 
MFs of (and best fits to) different pairs for stellar masses be- 
low 1 .OMq, both for the full STIS field and for areas limited by 
radial range. The results have been corrected for sample incom- 
pleteness and background contamination. In Liu et al. (2009) we 
combined the NGC 1818 MFs from WFPC2 and STIS. Using 
WFPC2, we detected more stars than on the basis of our STIS 
data for some mass ranges, so the slope of the MF in Liu et al. 
(2009), r - 0.46 + 0.10, is steeper than the slope obtained here 
(F - 0.28 + 0.04). However, they are both in agreement with that 
of the Kroupa (2001) IMF, F = 0.3±0.5, within the uncertainties. 

The MFs of NGC 1805, NGC 1818, and NGC 1868 show 
the same trend, as do the MFs of NGC 1831, NGC 2209, and 
Hodge 14. Figure |6] shows the change in trend most clearly, for 
our intermediate-age cluster pair Because NGC 1868 is much 
more compact than NGC 1831, its dynamical age is much older. 
This implies that it will be much more evolved dynamically than 
NGC 1831, hence resulting in a turnover at much lower masses 
(outside of our observational range). Similarly, NGC 2209 and 
Hodge 14 are an order of magnitude older (from both a stellar 
evolution and a dynamical point of view) than the Pair II clusters, 
hence exhibiting declining MFs in the mass range of interest here 
(see Fig.|2ll. 

We obtained all cluster MFs based on only single stars, ne- 
glecting unresolved binary and other multiple stars. As discussed 
in Liu et al. (2009), Kerber et al. (2006) analysed the effect of bi- 
narity on cluster MFs. They found that MFs with binary fractions 
of unity and 0% are identical within the observational uncertain- 
ties (cf. Liu et al. 2009 for NGC 1818), so that we can justifiably 
ignore the effects of binarity in the context of the low-mass sys- 
tem MFs derived in this paper (and in Liu et al. 2009). 
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Fig. 5. Low-mass cluster MFs for Pair I, NGC 1805 and NGC 
1818. 



4. Discussion 

Owing to both the large distances involved and observational 
limitations, it is difficult to obtain deep stellar MFs in ex- 
tragalactic environments. Kroupa (2001) studied the solar- 
neighbourhood IMF down to ~ O.OIM© and reported his often- 
quoted broken power-law distribution. Paresce et al. (2000) anal- 
ysed the MFs of a dozen Galactic GCs down to Q.IMq and de- 
rived a lognormal distribution below I.OMq. Much effort has 
also focussed on studying stellar MFs in the LMC (e.g.. Will et 
al. 1995; de Grijs et al. 2002b; Gouliermis et al. 2006a,b; Kerber 
et al. 2006; Da Rio et al. 2009). Recently, the LH 95 IMF ob- 
tained by Da Rio et al. (2009) reached down to 0.3 IMq, i.e.. 
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1868. 



14. 



much deeper than achieved previously beyond the Milky Way. 
Chiosi et al. (2007) studied the MFs of young SMC clusters 
down to O.IMq. In Liu et al. (2009), we probed - for the first 
time - the stellar MF in an extragalactic, low-metallicity envi- 
ronment down to O.15M0. In this paper, we applied the method 
of Liu et al. (2009) to all of our sample clusters, armed at assess- 
ing the evolution (if any) of the low-mass MF. 

Santiago et al. (2001), de Grijs et al. (2002c), and Kerber et 
al. (2006) studied the stellar MFs above I.OMq at different radii 
in these same clusters. They found that the slopes at different 
radii were significantly different, because significant degrees of 
mass segregation affect the PDMF shape, de Grijs et al. (2002b) 



studied the NGC 1805 and NGC 1818 MFs above LOM©. They 
found that the cluster MFs followed the Salpeter (1955) IMF 
quite closely, de Grijs et al. (2002c) studied the LFs of all sample 
clusters and concluded that the PDMFs of the clusters in each 
pair must be very similar. 

The relaxation time at the half-mass radius of compact star 
clusters can be written as (Meylan 1987) 



/r,h = (8.92 X 10^)- 



M, 



1/2 



<m) log(0.4M,ot/(m)) 



(1) 



where is the half-mass radius (in pc), Mtot the total cluster 
mass (in Mq), and (m) the typical mass of a cluster star (in Mg). 
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Table 2. Mass functions. 
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Fig. 8. Full MFs of the sample clusters arranged by pair: all open 
points are from this paper, based on STIS observations, while the 
solid points are from de Grijs et al. (2002c), based on WFPC2 
observations. The dashed line represents a broken power-law 
IMF, adopting the Kroupa (2001) high- and low-mass slopes, 
and the solid curve shows the original Chabrier (2003) lognor- 
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The timescale on which a cluster will have lost all traces of its 
initial conditions is well represented by its half-mass relaxation 
time. The dynamical properties of the Pair I clusters were dis- 
cussed by de Grijs et al. (2002b). They (see also de Grijs et al. 
2003) computed the half-mass relaxation time as a function of 
mass for NGC 1805 and NGC 1818. For stellar- masses below 
l.OMo these were > 300 and > 700 Myr, respectively. This is 
much longer than the clusters' ages, which implies that dynami- 
cal cluster evolution will not (or negligibly) affect the MFs below 
I.OMq at the clusters' half-mass radii. 

Based on the structural parameters of all clusters in Table 
[T] and their half-mass radii based on the WFPC2 and STIS ob- 
servations, we calculated the relaxation times of the other clus- 
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ters in pairs II and III for stellar masses below I.OMq. These 
timescales are > 2.0, 0.4, 3.9, and 0.52 Gyr for NGC 1831, NGC 
1868, NGC 2209, and Hodge 14, respectively, although dynami- 
cal evolution in the cluster core may proceed much faster Elson 
et al. (1987) calculated the relaxation time in the core and at 
the half-mass radius of NGC 1818 as log(f,-(0)/yr) = 8.2 - 8.8 
and log(fr.h/yr) - 9.0 - 9.7 (where the age range signifies the 
uncertainties due to the uncertain mass-to-light ratio); for NGC 
1831 the equivalent timescales were found to be log(fi(0)/yr) - 
9.1 - 9.5 and log(fr_h/yi") - 9.6 - 10.0. Compared to the clusters' 
ages (see Table 1), this implies that both clusters have undergone 
little to no significant dynamical evolution overall. 

From Figs.|5]|6l and|2]we conclude that the MP slopes for the 
entire STIS field of view and for smaller areas at different radii 
are identical within the uncertainties for a given cluster. There 
may be two reasons for this behaviour: (i) dynamical cluster evo- 
lution does not affect the MPs severely below 1 .QMq beyond the 
core region, although the relaxation times are shorter than their 
ages, at least for some clusters, or (ii) the STIS field is fairly 
small, which prevents us from detecting any differences in the 
MPs below I.OMq beyond the crowded centres. 

Because of the uniformity of the MPs derived from the full 
STIS field and from smaller areas at different radii, in Pig. [8] 
we compare the MPs of the clusters in each pair using the full 
STIS fields onlyQ In Liu et al. (2009), we adopted both bro- 
ken power-law and lognormal distributions to fit the NGC 1818 
MPs. Recent work supports this method (e.g.. Covey et al. 2008; 
Oliveira et al. 2009). In this paper, however, it has become clear 
that not all cluster MPs below 1 .OMq show an obvious turnover. 
Therefore, we adopted power laws to fit all cluster MPs below 
I.OMq. We include the relevant parameters, including those for 
the standard Kroupa (2001) IMP, in Table [3] The MP of NGC 
1805 shows the same result as that of NGC 1818 in Liu et al. 
(2009), i.e., their slopes are identical to the Kroupa (2001) IMP 
within the uncertainties, although the applicable mass ranges ex- 
tend to higher masses than the relevant Kroupa (2001) slope 
(this is likely due to statistical fluctuations; cf. Liu et al. 2009). 
NGC 1868 exhibits a broken power-law distribution similar to 
the Kroupa (2001) IMP and the complete IMPs of NGC 1805 
and NGC 1818 (de Grijs et al. 2002b; Liu et al. 2009). The MP 
of NGC 1831 does not show a turnover, but its slope in the mass 
range from 0.54 to I.I5M0 is identical to that of the Kroupa 
(2001) IMP for masses between 0.5 and l.OM©, again within the 
uncertainties. The MP slopes of NGC 2209 and Hodge 14 are 
also identical to those of the Kroupa (2001) IMP below I.OMq. 

In this paper, we extended our previous work (Liu et al. 
2009) to all LMC sample clusters. The MPs of the clusters in 
each pair exhibit identical slopes, and they are all also identical 
to the standard Kroupa (2001) IMP below l.OM©, independent 
of metallicity, particularly for the clusters in pair I (de Grijs et 
al. 2002b; Liu et al. 2009) and NGC 1868 in pair II (this pa- 
per), de Grijs et al. (2002c) studied the LPs of all sample clusters 
above 1 .OMq. We converted the LPs to MPs based on the Padova 
isochrones for the appropriate metallicities and combine both of 
our results to obtain complete PDMPs for all sample clusters 
(see Pig. [8] and Table |2]l. All MPs are similar to the standard 
Kroupa (2001) IMP, and the MPs of NGC 1831, NGC 2209, and 
Hodge 14 also match the Chabrier (2003) lognormal distribu- 



Table 3. Mass function slopes. 



Sample 


Low mass 


High 


mass 




olope (1 ) 


iviass {Mq) 


Slope (r) 


Mass (Mq) 


Kroupa 


0.3±0.5 


0.08-0.5 


— 1. j±U.J 


n ^ in 
U.J— l.U 


(2001) 






-1.3±0.7 


> 1.0 


NGC 1805 


0.21+0.10 


0.14-0.79 






NGC 1818 


0.21+0.02 


0.15-0.63 






NGC 1831 






-1.43±0.11 


0.54-1.15 


NGC 1868 


0.28±0.08 


0.25-0.63 


-1.42±0.29 


0.63-0.79 


NGC 2209 






-1.23+0.04 


0.23-0.60 


Hodge 14 






-1.24±0.16 


0.46-0.75 



^ Note that the gap in stellar masses for NGC 2209 between the STIS 
and WFPC2-based MFs is caused by incompleteness differences be- 
tween both sets of observations, and does not necessarily reflect a true 
lack of intermediate-mass stars in the cluster. 



tion, at least qualitatively although not in detail (as shown by the 
solid lines in the individual panels). 



5. Summary and conclusions 

We extended our pilot study in Liu et al. (2009) and used deep 
HST/STIS photometry of a carefully selected sample of rich, 
compact clusters in the LMC to derive their stellar MPs for 
masses below 1 .OMq, which we combined with the MPs of de 
Grijs et al. (2002c) above I.OMq to obtain complete MPs for 
all sample clusters. To our knowledge, together with Liu et al. 
(2009) this is the first time that anyone has probed stellar (clus- 
ter) MPs to this depth in an extragalactic, low-metallicity envi- 
ronment. Based on our STIS observations, the MPs of our com- 
pact clusters are all identical to the standard Kroupa (2001) IMP 
below I.OMq, within the uncertainties. 

The observations were taken beyond the crowded cluster 
cores and only included stars at the half-mass radii. In addi- 
tion, the relaxation timescales of low-mass stars are much longer 
than the equivalent periods for high-mass stars, so that dynami- 
cal evolution will not have affected the younger clusters' stellar 
mass distributions below I.OMq noticeably. We have therefore 
provided unprecedented insights into the IMP in a low-density 
and low-metallicity extragalactic environment. 
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